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SUMMARY 
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Basic  studies  on  the  effect  of  inhomogeneities  -  via  their  geometry,  size,  material  properties, 
distribution  ,  multiple  scattering  processes  and  statistical  correlations  in  position  and  the  effect  of 
aligned  versus  randomly  oriented  non-spherical  scatterers  has  continued  this  year.  Emphasis  has 
been  on  studies  of  the  second  moment  of  the  field  fluctuation,  namely  the  incoherent  intensity,  the 
effects  of  completely  correlated  (identical)  scatterers  versus  partially  correlated  scatterers  ( distinct ) 
has  been  studied  as  a  function  of  observation  angle,  non-dimensional  wavenumber,  and 
concentration.  Due  to  the  specific  form  of  the  expressions  for  the  incoherent  intensity,  the  effect  of 
correlations  is  to  decrease  the  value  of  the  incoherent  intensity.  Since  the  correlation  contribution  is 
also  very  much  angle  dependent,  the  effects  are  more  in  some  directions  than  in  others.  In  every 
case  both  vertical  and  horizontal  polarizations  are  studied  and  the  differences  in  the  results  are 
discussed.  The  multipole  amplitudes  of  the  effective  field  exciting  each  scatterer  are  explicitily 
calculated  leading  to  more  accurate  values  of  the  second  moment.  It  may  be  recalled  that  the 
effective  field  and  the  effective  propagation  consatant  are  calculated  in  the  Quasi  Crystalline 
approximation  and  include  the  effects  of  pair  correlation  and  sequential  multiple  scattering  terms. 
The  so  called  back  scattering  enhancement  is  also  investigated,  relative  to  the  experimental 
observations  of  Kuga  and  Ishimaru  and  the  theoretical  explanations  of  Tsang  and  Ishimaru  using 
'  cyclical '  diagrams. 
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INTRODUCTION 

In  this  paper,  a  systematic  study  is  made  of  first  order  contributions  to  the  second  moment 
or  incoherent  intensity  of  the  electromagnetic  field  propagating  in  a  discrete  random  medium.  The 
second  moment,  which  is  traditionally  defined  as  the  correlation  function  of  the  component  of  the 
field  fluctuations  in  any  direction  u,  denoted  by  ly  can  be  written  as 

Iu=  <u»(E-<E>)(E-<E>)*‘u>  =<(u*AE)(u‘AE)*> 
where  AE  =  E-  <E>  is  the  fluctuation  of  the  electric  field.  Since  the  field  fluctuations  can  be 
expanded  in  a  multiple  scattering  series  each  term  of  which  contains  sums  on  all  possible 
scatterers,  it  is  evident  that  we  can  divide  the  resulting  terms  into  two  sets,  one  of  which  involves 
considering  only  those  terms  in  which  the  same  scatterer  contributes  to  a  particular  order  term  in 
each  field  fluctuation  and  the  other  which  involves  distinct  scatterers.  This  latter  set  of  terms  will 
contribute  to  the  incoherent  intensity  only  if  statistical  correlations  between  scatterers  are  taken 
into  account.  The  first  category  of  terms  are  equivalent  in  spirit  to  radiative  transfer  theory  since  it 
is  essentially  the  intensity  of  the  field  scattered  by  each  scatterer  that  propagates  from  one  site  to 
another.  Even  for  this  set  of  terms,  positional  correlations  between  scatterers  should  be  taken  into 
account  at  volume  fractions  exceeding  a  few  %,  but  these  terms  contribute  to  the  incoherent 
intensity  even  if  correlations  are  neglected. 

Neglecting  correlations,  and  modeling  the  scatterers  as  dipoles,  Tsang  and  Ishimaru  have 
attempted  to  explain  the  back  scattering  enhancement  experimentally  observed  by  Kuga  and 
Ishimaru.  It  must  be  pointed  out  that  the  experiments  were  performed  at  wavelengths  comparable 
to  and  smaller  than  the  scatterer  size,  invalidating  the  dipole  approximation  used  in  the  theory, 
and  further  at  the  concentrations  considered  in  the  experiments,  correlation  effects  cannot  be 
neglected.  This  will  be  discussed  in  greater  detail  later. 

The  aim  of  the  study  presented  here  is  to  make  very  accurate  calculations  of  the  first  order 
contributions  to  the  incoherent  intensity.  By  accurate,  we  mean  that  the  calculations  will  be 
numerically  exact  up  to  first  order.  This  is  achieved  by  modeling  the  single  scatterer  response  by 


considering  the  full  T-  matrix  ( multipole  field ),  include  the  effects  of  pair  correlations  to  the  first 
order  term  of  the  incoherent  intensity,  and  lastly  the  field  incident  on  a  scatterer  is  modeled  as  an 
effective  exciting  field,  which  propagates  at  the  effective,  complex,  frequency  dependent 
wavenumber  K,  which  is  first  computed  in  the  Quasi  Crystalline  Approximation  (  QCA).  The 
multipole  amplitudes  of  the  effective  exciting  field  are  calculated  using  the  extinction  theorem. 
This  is  in  contrast  to  previous  calculations,  when  all  amplitudes  were  set  to  unity.  By  necessity 
the  calculated  intensity  is  normalized  by  the  incident  intensity,  the  illuminated  volume,  and  the 
distance  of  propagation  into  the  scattering  volume,  to  take  account  of  the  attenuation  that  the  wave 
suffers  before  it  is  incident  on  the  last  scatterer. 

Numerical  results  are  presented  in  various  forms,  as  a  function  of  frequency,  observation 
angle,  concentration,  polarization  of  the  scattered  field,  i.e.  vertical  and  horizontal  polarizations. 
An  effort  was  made  to  separate  the  contributions  to  first  order  scattering,  so  that  the  effect  of  pair 
correlations  is  clearly  visible.  To  first  order,  neglecting  correlations,  results  in  the  single 
scattering  approximation  to  the  incoherent  intensity.  The  effect  of  correlations,  is  to  actually 
reduce  the  intensity,  especially  close  to  the  forward  direction. 
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Incoherent  Scattered  Intensity 

The  total  scattered  intensity  is  directly  proportional  to  the  second  moment  of  the  scattered 
field.  It  is  known  that  the  total  scattered  field  is  a  combination  of  the  average  scattered  field  and 

t 

the  fluctuation  of  the  field  due  to  the  random  distribution  of  scatterers,  i.e.,  u  =  <u>  +  u,  the 
incoherent  component  of  the  scattered  intensity  can  be  obtained  as 

<uu  >  =  <  uu  >  -  <u><u>  (1) 

where  the  angular  brackets  denote  an  ensemble  average.  To  first  order,  that  is  taking  only  the 
single  scattering  contribution  to  each  scattered  field,  we  obtain 
<u  u  *>  =  <|u'l^>=  <Zuj*  Xujj>  -  X<Uj>*Z<Ulj> 

=  X  X  ^  <  I  Uj  I  -  XX  <up*<uit>  (2) 

j+k  k  j  j  k 
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In  the  above  equation,  the  subscripts  j  and  k  denote  the  scattered  field  from  the  j-th  and  k-th 
scatterers,  respectively.  Subscripted  angular  brackets  denote  a  conditional  ensemble  average  in 
the  usual  manner,  where  the  averaging  is  done  keeping  the  position  of  the  subscripted  scatterer(s) 
fixed.  For  scatterers  randomly  distributed  in  space,  the  ensemble  averages  in  Eq.  (2)  can  be 
written  by  integrating  over  the  random  positions  rj,  rj^,  etc.  of  the  scatterers.  Thus 
<  I  u  1 /  /  [(N-l/N)  <u*jUj{.>j}j  -  <Uj>j*  <Uk>kl 
+  Uq  1  <  1  uj  1 2>j  drj ,  (3) 

where  n^  =  N/V,  is  the  number  density  such  that  the  total  number  of  scatterers,  N->  «>,  the 
volume  in  which  they  are  distributed,  V— ><»,  but  n^  is  finite.  The  second  term  in  the  above 
equation,  which  is  proportional  to  n^  is  the  ordinary  single  scattering  approximation  to  the 
incoherent  intensity  and  the  magnitude  of  the  incoherent  intensity  in  any  direction  is  proportional 
to  the  bistatic  cross  section  of  a  single  scatterer.  The  first  term  proportional  to  n^^  is  due  to  the 
effect  of  positional  correlations  between  pairs  of  scatterers. 

For  sparse  concentrations,  the  scatterers  are  uncorrelated,  and 

<“j*“k>jk  =  <“j>j*<“k>k 
and  the  first  term  of  Eq.  (3)  vanishes. 

In  order  to  obtain  the  incoherent  intensity  for  higher  concentrations,  the  pair  correlation 
function  g(x)  is  taken  into  consideration.  Recently  Twersky  (1982)  has  modified  Eq.  (3)  for  a 
dense  distribution  of  scatterers  in  which.  The  incoherent  intensity  has  the  form: 

<lu'l2>  =  ng^  j|  fg(rj-rjj)  -  1)  <Uj>*j  <ujf>ij  drj  dr^. 

+  n^j  J<  1  luj  |2>j  dfj  (5) 

Equation  (5)  is  the  final  form  used  in  our  computation.  However,  a  knowledge  of  the 
coherent  field  is  still  required  since  the  average  scattered  field  <Uj>j  is  involved  in  the  formalism. 
In  the  usual  manner,  the  field  scattered  by  the  j-th  scatterer  is  expanded  in  vector  spherical 
functions,  with  coefficients  f^J,  which  arc  related  to  the  expansion  coefficients  of  the  exciting 
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field  via  the  T-  matrix.  Thus, 


and  denoting  cXjjJ  ,  the  coefficients  of  the  exciting  field 

f  j  =  T.,ri'Ot«*^ 

‘n  ^  nn  n 

The  average  scattered  field  <Uj>j  holding  the  j-th  scatterer  fixed  can  thus  be  expressed  as 

<Uj>j  =  Tnn-  <otnJ>j 

In  writing  Eq.  (5),  we  have  used  the  quasi-crystalline  approximation,  as  explained  in  the 
previous  section. 


(6) 

(7) 

(8) 


The  exciting  field  coefficients  ot^J  are  initially  unknown  in  Eq.  (5).  However,  the  average 
field  <otnj>j  exciting  the  j-th  scatterer  is  known  after  defining  an  effective  propagation  constant  K 
which  is  a  complex  (  K=Kj  +  i  K2).  Following  this  definition,  the  average  or  effective  exciting 
field  <cXj,J>j  can  be  written  as 

«Xnj>j  =  An  e *‘o*'‘j  (9) 

where  is  the  propagation  direction  of  the  incident  wave.  The  effective  propagation  constant  K 
is  obtained  by  solving  the  dispersion  equation  which  is  c  .plained  in  detail  in  our  previous  papers. 
The  unknown  effective  exciting  amplitude  An,  however,  can  be  solved  by  invoking  the  extinction 
theorem  and  is  generally  in  the  following  form  for  spherical  scatterers,  for  which  the  T-matrix  is 
diagonal 


An  =  [1  -3/2  c  (-1)"  I.Tn'An*]  [(ka’  )  (  Kfic  +  1)]  (10) 

provided  a  spherical  scatterer  in  an  effective  medium  is  considered.  In  Eq.  (8),  c  is  the  volume 
concentration,  a  is  the  radius  of  the  scatterer  and  k  is  the  wave  number  in  the  host  medium. 

If  we  substitute  Eqs.  (8)  and  (9)  in  (5),  we  obtain 
<  I  u  I  1 J  Ig(rj-r|;)  <«„j>j  V'nn"* 

.  (11) 

We  notice  that  in  Eq.  (11),  the  multiplication  of  the  T-matrix  and  the  effective  exciting  field  is 
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independent  of  the  integral  and  the  pair  correlation  function.  In  fact,  the  integral  of  the  pair 
correlation  function  turns  out  to  be  the  Fourier  transform  of  the  radial  distribution  function 
[g(x)-l].  This  fact  further  simplifies  the  numerical  computation  of  the  incoherent  intensity.  In 
order  to  investigate  the  contribution  of  the  incoherent  intensity,  we  calculate  two  major 
normalized  quantities  which  are  defined  as  follows: 

(i)  Normalized  incoherent  Intensity  <I  >  in  the  far  field 
<l'>  =  (VIo)(kr)2(v/V) 

=  2c(I  Tn-An-Yn-X  I  VA^.. ¥„..)*  F(ko,K,r)/  (K2/k)(ka)3(Z/a)  (12) 

where  are  the  normalized  spherical  harmonics,  (Yjjj^(j(0,(t))),  Z  is  the  distance  of  penetration 
of  the  incident  wave,  n'  is  the  index  representing  1  and  m;  r  is  the  radial  distance  from  the  center 
of  the  scattering  region  to  the  observation  point,  v  is  the  volume  of  a  single  scatterer  and  V  is  the 
whole  scattering  volume  (region).  The  function  F  in  Eq.  (12)  is  given  below 

F=  l+n^jJ  [r  g(x)-l]  exp  [  i(KkQ-k  r /r  )  •  X  ]  dx  (13) 


RESULTS  AND  DISCUSSION: 

fm  In  order  to  study  the  characteristics  of  the  incoherently  scattered  wave  intensity,  we  choose 

’  ■  electromagnetic  waves  as  probing  waves  simply  because  there  are  a  number  of  applications  in 

remote  sensing.  Recently,  there  have  been  efforts  to  explain  the  so  called  backscattering 

si* 

enhancement  both  experimentally  and  analytically.  Although  several  attempts  have  been  tried  to 
explain  the  phenomenon,  we  think  the  major  contribution  toward  a  locally  high  incoherent 
intensity  in  the  backscattering  direction  may  depend  more  on  Mie  scattering  and  less  on  the 
1 '  concentration  of  scatterers.  Besides,  the  properties  of  different  scatterers  do  affect  the  scattering 

\  characteristics  and  the  backscattering  enhancement  relies  heavily  on  the  albedo  of  a  single 

A* 

scatterer.  We  are  going  to  explain  this  later  by  using  several  of  our  numerical  results. 

By  sending  plane  electromagnetic  waves  through  scatterers  we  intend  to  find,  first,  the 
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angular  dependence  of  the  incoherent  intensity  and  the  influence  from  different  polarizations, 
e.g.,  vertical  and  horizontal  polarizations.  If  it  is  not  specifically  mentioned,  scatterers  are 
assumed  to  be  spherical  ice  particles  with  a  relative  dielectric  constant  e  =  3.168  embedded  in  air. 
Figure  1  presents  the  normalized  incoherent  intensity  versus  observation  (scattering)  angles.  The 
forward  scattering  angle  is,  in  our  case  0  °  and  therefore  the  backscattering  direction  is  at  180°. 
The  nondimensional  frequency  considered  is  0.6  which  is  equivalent  to  a  physical  frequency  of 
about  14  GHz  if  a  2  mm  particle  is  considered. 

After  the  electromagnetic  waves  travel  about  1  km,  we  can  see  that  the  scattered  incoherent 
intensities  in  all  directions  are  very  small  compared  to  that  of  the  coherent  wave  (or  the  coherently 
transmitted  wave).  As  a  matter  of  fact,  this  is  true  for  all  incoherently  scattered  intensities.  From 
the  field  investigation  as  well  as  from  the  controlled  experiments,  incoherent  signals  are  usually 
25-50  dB  or  even  more  down  from  the  coherent  signal.  But  as  we  discussed  in  the  previous 
section,  this  does  not  imply  that  the  incoherent  signals  are  unmeasurable. 

Taking  a  further  look  at  Fig.  2,  we  can  conclude  that  the  vertical  polarization  gives  more 
angular  dependence  of  the  incoherent  intensity  than  the  horizontal  one.  There  is  an  extremely  low 
intensity  (i.e.,  a  deep  minimum  in  the  curve)  that  occurs  at  90°  at  a  ka=0.6.  This  phenomenon 
happens  again,  however,  at  a  higher  observation  angle  of  125°  when  the  frequency  ka  is  raised 
to  2.  There  is  no  polarization  difference  at  the  forward  and  backscattering  directions  for  the 
incoherent  intensities. 

Numerical  results  for  the  normalized  incoherent  intensity  versus  frequency  for  different 
concentrations  at  the  backscattering  direction  are  presented  in  Fig.  3.  The  scattered  intensities  are 
vertically  polarized.  As  can  be  seen  from  the  figure,  locally  low  intensities  happen,  for  all 
concentrations,  at  ka=1.5  which  is  a  resonant  frequency  for  ice  particles.  Besides  the  magnitudes 
of  the  intensities,  one  is  able  to  find  that  the  minimum  is  wider  when  the  concentration  is  smaller. 
The  difference  in  the  magnitudes  of  the  intensities  is  quite  noticeable  in  the  frequency  range  from 
ka=0.7  to  1.5. 
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In  Figs.  4  and  5,  we  compare  the  backscattered  intensity  calculations  with  and  without  the 
effect  of  pair  correlations.  These  results  tell  us  that  if  the  intensity  is  calculated  without 
considering  the  pair  correlation  function  when  the  concentration  becomes  even  moderately  high, 
i.e.,  5%,  one  is  able  to  see  the  difference  in  the  magnitude  particularly  in  the  low  frequency 
range.  This  is  valid  for  both  vertical  and  horizontal  polarizations. 

The  scatterers  in  the  host  medium  are  actually  excited  by  waves  with  the  effective  exciting 
amplitudes  which  result  from  different  orders  of  scattering.  One  way  to  consider  the  order  of 
scattering  is  by  perturbing  the  incident  wave  amplitude  which,  in  our  case,  is  unity.  Therefore, 
in  Fig.  6,  we  present  results  from  calculating  incoherent  intensities  using  different  exciting 
amplitudes  and  the  difference  is  quite  noticeable  when  the  frequency  is  high,  especially  at  the 
resonant  scattering  frequency. 

Figure  7  is  another  representation  of  the  backscattered  intensity  versus  frequency. 
However,  the  dependence  of  the  intensities  on  concentrations  is  clearer  using  the  backscattering 
coefficient  computation.  The  bistatic  scattering  coefficient  has  a  strong  angular  dependence  under 
different  frequencies  and  is  shown  in  Fig.  8. 

In  order  to  understand  the  effects  of  multiple  scattering  on  the  incoherent  intensity  and  its 
final  characteristics,  the  study  of  the  scattering  response  from  a  single  scatterer  proves  to  be  very 
informative.  First,  we  examine  the  differential  scattering  cross  section  for  an  ice  sphere  at  low 
and  high  frequencies.  The  extremely  low  intensity  occurring  at  90°  in  Fig.  1  and  can  be 
explained  by  the  polar  plot  of  bistatic  scattering  cross  sections  in  Fig.  9.  At  high  frequencies,  the 
intensity  has  a  strong  lobe  in  the  forward  direction  and  is  highly  directional  in  the  range  of  about 
±10°.  However,  if  a  perfect  conducting  sphere  is  considered  as  a  scatterer,  at  low  frequencies, 
we  can  see  from  Fig.  1 1  that  the  forward  scattering  is  about  four  times  smaller  than  the 
backscattering.  Furthermore,  the  scattering  cross  section  is  not  zero  at  90°.  This  means  that  the 
material  properties  of  different  scatterers  dominate  the  final  scattering  characteristics.  Therefore, 
even  at  the  same  observation  angle,  a  particular  phenomenon  that  occurs  for  one  group  of 


scatterers  may  not  occur  for  another  group  of  scatterers  with  different  properties. 

The  characteristics  of  the  multiply  scattered  incoherent  intensity  has  its  similarity  with  that  of 
the  differential  cross  section  of  a  single  scatterer  and  can  be  again  found  in  Figs.  12  and  13  when 
compared  with  Figs.  1  and  2.  Figure  14  is  a  different  representation  of  the  polar  plot  in  Fig.  10 
but  detailed  magnitudes  with  respect  to  the  observation  angles  are  given.  In  addition,  the 
difference  between  the  vertical  and  horizontal  polarizations  is  shown  in  the  figures. 

Different  scatterers  have  a  different  albedo  when  illuminated  by  waves  with  different 
frequencies.  Since  the  incoherent  scattering  characteristics  depends  heavily  on  the  albedo  of  a 
single  scatterer,  Figs.  15  and  16  present  the  differential  scattering  cross  sections  versus  scattering 
angles  for  an  ice  particle  in  air,  a  perfect  conducting  sphere  in  free  space  and  a  glass  sphere  in 
distilled  water  all  excited  by  plane  waves  of  two  different  frequencies,  i.e.,  ka=0.6  and  7.28. 
Indeed,  we  can  see  the  property  dependence  and  some  behavior  related  to  the  incoherent 
scattering  may  be  explained  from  these  figures. 

The  differential  cross  section  of  a  single  scatterer  also  has  a  noticeable  frequency 
dependence.  In  Fig.  17,  we  observe  that  the  resonant  scattering  is  different  for  different 
scatterers.  For  an  ice  sphere  in  air  it  happens  at  ka=1.5  while  it  shifts  to  a  higher  ka  (=1.6)  for  a 
glass  sphere  in  water.  This  implies  that  extra  information  is  always  available  if  we  pay  attention 
to  the  details  of  the  frequency  spectrum  of  the  incoherent  intensity. 

Finally,  we  want  to  say  something  about  the  effect  of  the  pair  correlation  function.  In  order 
to  tell  the  importance  of  its  effect  on  the  final  scattered  intensity,  we  simply  calculate  the  function 
F  which  appears  in  Eq.  (12)  and  has  been  defined  in  Eq.(13).  As  can  be  seen  from  Eq.  (12),  it 
involves  a  Fourier  transform  of  the  pair  correlation  function  and  it  contains  an  effective 
propagation  constant  K  and  hence  depends  on  the  properties  of  the  scatterers,  the  concentration, 
frequency  and  angle  of  observation.  However,  in  Figs.  18,  19  and  20,  we  can  see  that  the 
Fourier  transform  of  the  pair  correlation  function  dominates  the  scattering  response  particularly  in 
the  low  frequency  range  and  in  the  forward  direction.  In  the  high  frequency  range,  it  does  not 


affect  the  scattered  intensity  ranch.  Also,  we  observed  that  the  intensity  decreases  after  a  volume 


fraction  c  =  15%  ,  which  is  also  a  fact  pointed  out  in  Kuga's  experiments.  One  is  therefore  able 


to  conclude  that  when  the  frequency  is  high  the  magnitude  of  the  scattered  intensity  is  more 


dominated  by  Mie  scattering  as  shown  in  Fig.  16.  At  low  frequency,  the  backscattering 


enhancement  cannot  be  observed  especially  for  low  concentrations  of  scatterers.  Even  if  the 


concentration  of  scatterers  is  increased,  no  backscattering  enhancement  at  the  low  frequency 


range  can  be  observed. 
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Normalized  incoherent  intensity  vs.  nondimensional  frequency 


Figure  6.  Normali-ed  incoherent  intensity  vs.  nondimens ional  frequency. 


Polar  plot  of  tistatic  scattering  cross  section. 
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Figure  17.  Bistatic  scattering  cross  section  vs.  nondimensional  frequen 
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